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ABSTRACT 

I Space telescopes such as EChO (Exoplanet Characterisation Observatory) and JWST 

(James Webb Space Telescope) will be important for the future study of extrasolar 
planet atmospheres. Both of these missions are capable of performing high sensitiv- 

I ity spectroscopic measurements at moderate resolutions in the visible and infrared, 

which will allow the characterisation of atmospheric properties using primary and sec- 
onda ry transit spectro scopy. We use the NEMESIS radiative transfer and retrieval 

I tool (|lrwin et al . 2008; Xee et al.ll2012f ) to explore the potential of the proposed EChO 

mission to solve the retrieval problem for a range of H2-He planets orbiting different 

I stars. We find that EChO should be capable of retrieving temperature structure to 

~200 K precision and detecting H2O, CO2 and CH4 from a single eclipse measurement 
for a hot Jupiter orbiting a Sun-like star and a hot Neptune orbiting an M star, also 

\ providing upper limits on CO and NH3. We provide a table of retrieval precisions for 

these quantities in each test case. We expect around 30 Jupiter-sized planets to be 
observable by EChO; hot Neptunes orbiting M dwarfs are rarer, but we anticipate 

\ observations of at least one similar planet. 

Key words: Methods: data analysis - planets and satellites: atmospheres - radiative 
transfer 



wavelength gives the emission spectrum of the planet's day- 
side, which as well as providing information about absorbing 
gases can place constraints on the temperature structure. 

Whilst missions like EChO and JWST would signifi- 
cantly advance our capability to perform this kind of mea- 
surement, the degree of improvement in atmospheric re- 
trievals is highly dependent on the spectral range, resolu- 
tion and signal to noise of the instrument. Solving the in- 
verse problem (jRodgerd |200Q| ) for planetary atmospheres 
is complex and is often hampered by insufficient infor- 
mati on in spe c tra to break degeneracies in the solutions. 
Tess envi et al.l (|2012l l explore the number of measurements 
required to observe primary and secondary transit spectra 
with a dedicated space mission for different planet cases, but 
they do not consider retrievability for different atmospheric 
scenarios. This paper explores the potential of an EChO-like 
telescope to solve the retrieval problem for a range of plan- 
ets orbiting different stars; the results are also applicable to 
any telescope with a similar spectral range (0.4 — 16 M^m), 
resolving power and noise level. 

We use the NEMESIS ra d iative transfer and retrieval 
tool developed by llrwin et al.l (| 20081 ) to generate synthetic 
spectra for model giant planets under different conditions, 
with the expected EChO spectral range and resolution. We 



. 1 INTRODUCTION 

The future launches of space telescopes such as EChO (Ex- 
I oplanet Characterisation Observatory, iTinettil I2OI2I ) and 
p JWST (James Webb Space Telescope) have the potential to 
introduce a new era in the study of extrasolar planet atmo- 
spheres. Designs for both missions incorporate the capability 
to perform high sensitivity spectroscopic measurements at 
moderate resolutions in the visible and infrared. Observa- 
tions will be made of exoplanets that transit their host star 
from the perspective of an observer on Earth; properties 
of their atmospheres can then be determined by comparing 
the in-transit and out-of-transit fluxes of the system over a 
range of wavelengths, a technique called transit s pectroscopy 
(ICoustenis et al.ill997l ; ISeager S asselovllioool ). This tech- 
nique translates the small additional reduction in flux caused 
by the planet's atmosphere as the planet crosses the stel- 
lar disc to a measurement of atmospheric opacity at each 
wavelength, which in turn provides information about the 
atmospheric scale height, aerosols and absorbing gases in 
the atmosphere. When the planet is eclipsed by the star, 
the difference between in and out of transit fluxes at each 
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then add noise to the synthetic spectra and feed them back 
into NEMESIS to retrieve the atmospheric properties of the 
model planet. Comparing the retrieved atmospheric state 
with the original model indicates whether EChO could pro- 
vide sufficient information for NEMESIS to correctly solve 
the retrieval problem for a real planet under similar condi- 
tions. 

NEMESIS, the Non-linear optimal Estimator for Multi- 
variateE spectral analysis, was originally developed to anal- 
yse data from the CIRS infrared instrument on the Saturn 
probe, Cassini. It has since been developed to work for any 
solar system body, and more recently to simulate primary 
transit and seco ndary eclipse spectra for extrasolar planets 
(|Lee et al.l l2Q12h . It is a particularly suitable tool for this 
task because of its proven versatility, and also its efficient 
approach to solving the radiative transfer equation. 

NEMESIS utilises the correl ated-k approximation 
(jGoodv Yundll989l : [Lacis Oinaslll99lh to calculate fast 
forward models, which speeds up the retrieval process. The 
correlated-k approximation allows absorption coefficients k 
to be ranked in order of strength across a spectral inter- 
val and pre-tabulated, and relies on the assumption that 
absorption line strengths are well- correlated between model 
atmospheric layers, i.e. lines that are strongest in the lowest 
atmospheric layer are also strong in the layer above. This 
approach reduces calculation time over the line-by-line ap- 
proach by reducing the number of ordinates over which the 
integration takes place, and the fast forward model calcula- 
tion i s coupled with an efficient optimal estimation scheme 
re.g.. lRodge73l20QOl ). The user provides NEMESIS with an 
initial guess and an associated error, the a priori solution, 
which acts to prevent overfitting and stops retrieval solutions 
from becoming unphysical. The best-fit solution is found by 
an iterative approach in which a forward model and the 
radiance derivatives with respect to the parameters to be 
retrieved are calculated at each step, and compared with 
the input spectrum; the next spectrum can then be calcu- 
lated based on the radiance derivative and the value of a cost 
function, which represents the difference between the mea- 
sured and synthetic spectra together with the deviation from 
the a priori solution. The optimal solution is achieved when 
the cost function is minimised. For further details about the 
structure of NEMESIS an d its use for retrievals of extraso- 
lar pla net atmospheres, see llrwin et al.1 (|2QQ8l ) andlLee et al.l 
(|20l2). 



2 ECHO 

At the time of writing, the proposed structure of the EChO 
telescope (as defined in the ECh O Science Requirements and 
Payload Definition Docu ments: lEChO Science Studv TeamI 
I2OIII : lEChO TeamI l201lh is a Cassegrain telescope with a 
flat folding mirror and a set of dichroic mirrors for spectral 
separation of the light. The effective collecting area is 1.131 
m^ and the field of view 20"x20". The light wih be split 
into six or seven channels, covering the wavelength range 
between 0.4 and 16 ^im. The visible channel detector, cov- 
ering the 0.4 — 0.8 ^im range, is likely to be a silicon-based 
CCD; the short to medium infrared channels, 0.8 — 1.5, 1.5 — 
2.5 and 2.5 — 5 ^im, will use HgCdTe detectors; the long- 
wave infrared channels, 5 — 8.5, 8.5 — 11 and 11 — 16 p.m, will 



Planet 


Mass (xl024 kg) 


Radius (km) 


Hot Jupiters 


1800 


75000 


Warm Jupiter 


1800 


75000 


Hot Neptune 


180 


30000 



Table 1. Masses and radii of the model planets in this study. 



probably be Si: As detectors, althou gjh other options are be- 
ing investigated (lEChO Teamll2Qllh . The target spectral re- 
solving power is R~300 at wavelengths less than 5 ^im, and 
R'^30 longwards of this. Possible descopes and trade-offs 
that might reduce this wavelength range and/or resolving 
power are not discussed here and are beyond the scope of 
this paper. 

There are two major advantages to EChO over JWST 
for observations of extrasolar planets; the first is that EChO 
is dedicated to exoplanet observations, meaning that more 
time is available for revisiting interesting targets and open- 
ing up the possibility to study temporal variation in exo- 
planet atmospheres (for targets around quiet stars); the sec- 
ond is that a single instrument covers the full wavelength 
range between 0.4 and 16 |J,m, avoiding the necessity of 
stitching together temporally dislocated results from differ- 
ent instruments to make a full spectrum. Stellar activity 
and instrument systematics make accurate stitching of spec- 
tra in different wavelength ranges taken at d ifferent times 
a com plicated, sometimes impossible, process: iTinetti et al.l 
(I2OIOI ) discuss this problem in the context of averaging mul- 
tiple observations taken within the same spectral range. 
JWST will launch earlier with a significantly larger mirror 
(an effective collecting area of 25 m^ versus EChO's 1.131 
m^) but its exoplanet science will be limited to a few select 
targets, whereas EChO hopes to characterise at least 100 
atmospheres. 



3 MODEL PLANETS AND SYNTHETIC 
SPECTRA 

In order to examine the level of information about exoplanet 
atmospheres that would be available from EChO, we use 
NEMESIS to generate a series of synthetic spectra based on 
a range of model planets. The parameter space that could 
be explored is very large, and we have reduced this to a few 
cases that we think best represent EChO's main targets. In 
this paper, we only discuss H2-He gas giant planets with 
trace amounts of H2O, CO2, CO, CH4 and NH3, of approxi- 
mately Jupiter and Neptune size. In future papers we intend 
to investigate effects of alkali metals, TiO/VO and haze on 
the visible spectra of gas giants, and extend our analysis 
to the super-Earth regime. All limits quoted f or EChO are 
taken from the Payload Definition Document (|EChO TeamI 
l201lh unless otherwise specified. 

We consider three different types of Jovian planet; a hot 
Jupiter-size planet, a hot Jupiter-size planet with a temper- 
ature inversion and a warm Jupiter-size planet all orbiting 
a Sun-like star at 35 pc. We use the Kurucz model solar 
spectrun{3, which has been smoothed to the EChO spectral 
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resolution, and assume a stellar radius of 695000 km. The 
value of 35 pc was chosen as it is just under half the stated 
distance to the faintest G-type target (75 pc), and is 15 pc 
greater than the distance to the brightest G-type target, so 
represents a dista nce which should be well within the capa- 
bilities of EChO. iBelu et alJ (|201lh state in their study for 
JWST that they expect ^^10 transiting hot Jupiters out to 
50 pc, so extrapolating to 75 pc we conclude that around 30 
of these objects should be within observable range. 

We also include a hot Neptune-size planet around an 
M-type star. At 35 pc, there is insufficient signal to observe 
a planet the size of Neptune around a sun-like star with a 
single transit. We place the M star at 6 pc, which in the 
baseline case for EChO is the greatest distance over which a 
transit around an M star is observable, with a goal of extend- 
ing this to 16 pc. Fewer studies are available of the expe cted 
occurrence rate for these planets; Bonfils et al.l (|201lh es- 
timate that 3 out of every 100 M dwarfs would have a hot 
Neptune in orbit, and there are around 1000 bright M dwarfs 
(JjlO) known within 16 pc (jLepine Gaidosll201lh . Despite 
the low occurrence rate, then, the large number of M dwarfs 
indicates that at least one hot Neptune should be observable 
by EChO, assuming a transit probability of at least 1/30. 
We use the Kurucz model spectrum for M5V, and assume 
a radius of 97995 km (equal to that of Proxima Centauri, 
an M5.5V star, as given in iDemorv et all l2009l ). Primary 
(transmission) and secondary (eclipse) transit spectra are 
obtained for all of these cases. 

We fix the temperature profile and bulk H2-He com- 
position of the planet's atmosphere, but allow the abun- 
dance of five trace gas species H2O, CO2, CO, CH4 and 
NH3 to vary independently of each other over the concen- 
tration range 0.1 — 1000 ppmv. The retrievability of different 
temperature profiles is implicitly tested by considering the 
different planet cases. The bulk at mosphere compos ition is 
assumed to be 90% H2 and 9.9% He: lLee et aP (I2OI2I ) found 
that variation in the He/H2 ratio had a negligible effect on 
the retrieval of temperature and molecular abundances for 
HD 189733b, which is confirmed in the present study. We 
compare retrieved temperature and gas abundances for hot 
Jupiter model spectra with 90% H2 and 9.9% He, retrieved 
firstly assuming the correct H2/He ratio and secondly as- 
suming 80% H2 and 19.9% He, and find negligible differ- 
ences. 

Very few assumptions have been made about the ex- 
pected levels of trace gases in each of these atmospheres, as 
the objective of this study is to investigate whether features 
due to these gases, if present, would be detectable in each 
case. Initially, we assume a constant volume mixing ratio 
(VMR) as a function of altitude for all gases, which is the 
simplest possible scenario. In Section [531 we take a brief look 
at a more 'realistic' vertical profile for CH4 and examine how 
our findings might be affected by changes in the VMR as a 
function of altitude. 

We generate between 100 and 500 synthetic transit and 
eclipse spectra for each test scenario, each with different 
uniformly distributed trace gas VMRs. We then retrieve the 
VMRs (and temperature as a function of atmospheric pres- 
sure in the eclipse case) from each synthetic. We calculate 

^ http:/ /kurucz. harvard.edu/stars/ 



Gas Source 

H2O HITEMP2010 (Rothmane^aL 2010) 

CO2 CDSD-1000 fTashkun et al. 2003) 

CO HITEMP1995 fRothman et al. 1995) 

CH4 STDS rWenger Champion 1998) 

NH3 HITRAN2008 fRothman. L. S. et al. 2009) 



Table 2. Sources of gas absorption line data. 

a reduced parameter for each retrieval, which represents 
the goodness of fit within the measurement error and is the 
divided by the number of degrees of freedom in the re- 
trieval. This value should be '-^l for a good retrieval. 

3.1 Gaseous absorption 

Line strength information for H2O, CO2, CO and CH4 
and colli sion-induc e d abs orption for H2-H2 and H2-He are 
all as in iLee et aT ] (I2OI2I ) and references therein. In addi- 
tion, we include NH3 wi th line strengths taken fro m the 
HITRAN2008 database (iRothman. L. S. et al.1 12OO9I ). The 
limited information available for some of these gases, par- 
ticularly our omission of high-temperature line informa- 
tion for NH3, may result in overly pessimistic predictions 
about retrievability of some constituents; in the future, 
we hope to exploit new infor mation from sources such as 
iTennvson Yurchenkol (|2012l ) to improve the completeness 
of gas absorption information in our radiative transfer mod- 
els, which will in turn increase the reliability of our pre- 
dictions. The sources of line data for each of the five trace 
absorbers are summarized in Tabled 

3.2 Spectral sensitivity 

The level of sensitivity to different atmospheric properties 
in transmission and eclipse (secondary transit) geometries 
is indicated by the derivative plots in Figure [T] It can be 
seen that there is equal sensitivity to temperature across a 
wide range of altitudes for the eclipse case when the whole 
spectrum is considered, but the sensitivity to temperature 
in transmission is heavily skewed towards the atmosphere 
above the 10-mbar level, because the atmosphere is opaque 
out to higher altitudes in transmission. This means that 
there is insufficient information to retrieve temperature as a 
function of pressure in transmission geometry. There is very 
little sensitivity to CO and NH3 in either geometry, but 
there are several spectral regions in which H2O, CO4 and 
CH4 have a measurable effect on the spectrum. The narrow 

^ The number of degrees of freedom is given by rimeasurements — 
^parameters — 1- When retrieving a continuous profile with some 
correlation length, it is difficult to exactly define the number of 
parameters. Here, we simply use the number of atmospheric levels 
for a continuous retrieval, plus the number of gas scaling factors 
retrieved. There are fewer parameters than this in reality, so our 
estimates of x^/n are conservatively large, but since the number 
of measurements (605) far outweighs the maximum number of 
parameters (55) we do not severely overestimate the x^, and it is 
clear that even with a slight overestimation we can achieve jn 
values of order 1. 
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Figure 1. Sensitivity of a hot Jupiter planet to various atmospheric parameters. The top plots show the derivative of the squared radius 
ratio with respect to a IK change in temperature at each altitude level (left) and with respect to a multiplication of the whole vertical 
profile of each gas (right). Similar plots are shown below for the change in flux ratio in eclipse. 



regions of sensitivity to CO and NH3, at 5 and 4 ]iin respec- 
tively, are heavily contaminated by absorption lines of the 
other gases, producing degeneracies in the solution. A key 
region of the spectrum for CO2 is the 15 ^im band; maximum 
sensitivity occurs at the centre of this band in transmission 
geometry, and in the wing of the band in eclipse. 

3.3 Noise Calculation 

We assume that the instrument is entirely photon-limited, 
and estimate the noise based on the square root of the num- 
ber of photons per channel during the transit. Whilst we 
may expect noise to be as much 2 — 3x photon- limited, this 
could easily be reduced to the photon limit by averaging 
over a greater number of transits, provided the noise is un- 
correlated. Star spots are likely to be problematic for targets 
orbiting M stars as these are often active, and changing con- 
figurations of spots could present difficulties when averaging 
over multiple transits. The presence of correlated noise and 
systematic noise related to the instrument or star would be 
a severe problem, but the effect is difficult to quantify for 
EChO at this stage of the design and is beyond the scope 
of this paper. The equations used to calculate the noise per 
channel are given below. 

nx = Ix7T {r,/D,f (X/hc) (X/R) QEr]t (1) 

where nx is the number of photons received for a given wave- 
length A, Ix is the spectral radiance of the stellar signal, 
is the stellar radius, is the distance to the star, h 
and c are the Planck constant and speed of light, R is the 
spectral resolving power, Aefr is the telescope effective area, 
QE is the detector quantum efficiency, 77 is the the through- 
put and t is the exposure time. The resolving power, effec- 
tive area, quantum efficiency (assumed 0.7 for all channels) 



and throughput per channel are all taken from the EChO 
Payload Definition Document (jEChO Teamll201lh . The esti- 
mated throughput 7] of each channel is given as 0.191 (0.4 — 
0.8 i^m), 0.284 (0.8—1.5 ^im), 0.278 (1.5—2.5 |^m), 0.378 
(2.5—5 ^im), 0.418 (5—8.5 |am), 0.418 (8.5—11 |j,m), 0.326 
(11 — 16 \im). The effective telescope area is 1.131 m^ The 
exposure time is estimated based on a typical transit dura- 
tion of 3 hours for the Sun-like star and 1 hour for the M 
dwarf, with an 80% duty cycle. 

To calculate the noise on the stellar signal, the square 
root of the number of photons is taken and then Equation [T] 
is reversed to obtain this in terms of spectral radiance, which 
gives the absolute noise on the observed radiance. We use 
the differential chain rule to calculate the noise in the ratio 
of in/out of transit fluxes, and assume that the planet signal 
is sufficiently small that /in^/out, where / is the radiance. 
The error in the flux ratio is then given by Equation [21 

erratic = ]j 2 (2) 

In the transmission case, the radius ratio 
(Rpianet/Ri.)'^ = 1 - /in/Zout- In the cclipse case, the 
planet-star flux ratio Fpianet/K^ — ^out/^in — 1. It can be 
seen from both of these equations that the error in the 
in/out of transit radiance ratio is also the error in the 
radius squared and flux ratios for the two transit cases. 

The error calculated in this way is taken to be the la 
value for Gaussian noise on the spectrum, so the final error 
added to each band in the synthetic spectrum is calculated 
by generating a series of Gaussian-distributed random num- 
bers with a mean of zero and a equal to the calculated noise. 

The scenario described above is of course a perfect sce- 
nario; whilst the aim for a telescope like EChO is to ap- 
proach the photon noise limit as closely as possible, the 
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Figure 2. An example of a fit to a noisy synthetic hot Jupiter 
ecHpse spectrum. The residual is shown in the lower panel. The 
l-cr error bar is shown in dark gray shading, and the 2-cr error 
in light gray. The joined crosses in the top panel show the noisy 
synthetic, and the solid black line shows the fitted spectrum. 



reality is likely to incl u de re d noise and instrument sys- 
tematics. iGibson et al.1 (|201lh compare the error bars on 
a HST/NICMOS transmission spectrum obtained by fully 
taking into account instrument systematics using Bayesian 
methods with work assuming white noise only, and find that 
the errors are inflated by a factor of ^3. This is an extremely 
conservative estimate, and since NICMOS systematics were 
particularly severe we may assume that the same increase 
in error for EChO represents a worst-case scenario. 



4 PLANET TEST CASES 
4.1 Hot Jupiter 

The first planet we consider is a hot Jupiter-size planet or- 
biting a sun-like star. This planet is similar to HD 189733b, 
which has been observed using the Hubble Space Telescope; 
evidence for H2O, CO2, CO and CH4 in the atm osphere 
of HD 189733b has been prese nted bv I Swain et al. 



iMadhusudhan SeageJ (|2009') and iLee et al.l (I2Q12V NH3 



has also been included as a possible constituent in this study 
as it has been detecte d in ^1000 K brown dwarf atmospheres 
(ISaumon et aLl lioOOl) which may have similarities with hot 
exoplanet atmospheres. To demonstrate sensitivity to dif- 
ferent vertical temperature profiles, we include two variants 
of the hot Jupiter case with and without a temperature in- 
version above the tropopause; the case with a temper ature 
inversion is similar to HD 209458b (iDesert et al.ll20Q9l ). Ini- 
tially, we assume that the atmosphere is well-mixed and the 
composition is constant as a function of altitude. 

The scale height of the atmosphere is calculated based 
on the mass and radius provided (see Table [1]) , from which 
we obtain the gravitational acceleration at the 'surface' of 
the planet. In the case of a gas planet, we assume that this 
'surface' lies at a pressure level of 10 bar, as the atmosphere 
is opaque at pressures greater than this. The atmospheric 
scale height for the hot Jupiter varies as a function of tem- 
perature, but is of order a few hundred km. 



The initial test incorporates 500 secondary eclipse spec- 
tra with different trace gas abundances. We utilise the whole 
of the anticipated EChO spectral range, and assume that 
measurements were made during a single transit. An exam- 
ple noisy synthetic spectrum, with best-fit model, is shown 
in Figure [2I 

We retrieve temperature as a function of pressure and 
altitude-independent volume mixing ratios for each of the 
five trace absorbers. The results are presented in Figures [3] — 
[4l The rainbow colours of each line/point corresponds to the 
reduced (x^ divided by the number of degrees of freedom 
n) of the fit in each case, with purple points having the 
lowest /n and red points the highest. All but a very few 
retrievals have a /n below 1.2. It can be seen that in 
general NEMESIS can successfully retrieve the temperature 
profile and VMRs of H2O, CO2 and CH4 in all cases; it is 
however not capable of retrieving CO and NH3 where the 
VMR of these gases falls below 10~^, so only a rather high 
upper limit would be achievable for these gases. The retrieval 
of H2O VMRs starts to become unreliable for VMRs below 
5x10"^. 

The quality of the temperature retrieval is good in the 
troposphere and stratosphere, but at higher altitudes and 
in the very deep atmosphere results are comparatively poor. 
This is due to degeneracies between the temperature re- 
trieval and the retrieval of H2O and CO2 abundances; when 
retrieving several quantities, there are sometimes multiple 
atmospheric scenarios that result in spectra with an equally 
good fit to the data, so the retrieval can converge on an 
incorrect solution. This is demonstrated in Figure [S] low 
VMRs of H2O and CO2 result in the retrieval of a higher 
mesospheric temperature than the 'real' temperature. In- 
formation about temperature at a particular altitude is ob- 
tained from fitting the shape of a H2O or CO2 absorption 
feature, so a lower abundance of either gas, corresponding 
to a smaller feature, may result in a reduction of available 
information from the feature. This could mean that the re- 
trieved temperature relaxes back towards the a priori^ which 
in this case is higher than the 'real' temperature. 

It can be seen in these plots that in cases where the 
retrieved value is far from the input the error bars can be 
small, e.g. for small VMRs of CO2. This is because NEME- 
SIS estimates the error on each retrieved parameter based 
only on the errors in the measured spectrum and in the a 
priori^ which does not account for degeneracy. It is there- 
fore crucial to explore the possibility of degenerate solutions 
by varying the a priori and performing the retrieval multi- 
ple times, which ensures that the interde pendence of solu - 
tions for different variables is understood. iLee et al ] (120121 ) 
do this, and they also examine the variation in retrieved 
temperature as a function of x^ /n for different gas VMRs. 
This provides a more robust estimate of the error in the 
retrieved temperature. 

As well as exploring degeneracies, it is also important 
to test the effect of different a priori assumptions on the re- 
trieved temperature profile. This test indicates the regions 
of the atmosphere for which information is available in the 
spectrum, as outside these regions the retrieved profile re- 
laxes back to the a priori. We test this by retrieving the first 
100 spectra from the previous test, but assuming three dif- 
ferent a priori temperature profiles. The results of this test 
are shown in Figure [6] It is clear that the standard deviation 
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Figure 3. Eclipse retrieval results for a hot Jupiter without a temperature inversion, colours correspond to different values, increasing 
from purple to red. The dashed black line on the temperature plot is the a priori temperature provided to NEMESIS; the solid black 
line is the input temperature profile. The black lines on the VMR plots are 1:1 correspondence lines. Some of the temperature retrievals 
with the very highest are poor, indicating that in these few cases NEMESIS has failed to adequately fit the spectrum. 




in retrieved temperatures decreases as the a 'priori becomes 
closer to the correct solution. The mean retrieved tempera- 
ture is correct for all three cases below the tropopause and is 
close to the correct value at pressures greater than 1 mbar; 
the mean retrieved temperature is correct everywhere for 
the case where the shape of the a priori profile is the same 
as the shape of the input profile. For the temperature at al- 



titudes above the 1-mbar level there is a lack of information 
available in the spectrum, so the retrieval coverges on the 
a priori value. From this test, we see that a better a pri- 
ori guess does slightly improve the accuracy of the retrieved 
temperature, but a comparatively poor a priori assumption 
can still produce a valid solution provided there is informa- 
tion in the spectrum. The jn is very similar for all a priori 
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Figure 5. Retrieval results for temperature; colours in the upper 
and lower plots correspond to different input VMRs of H2O and 
CO2 respectively, with low to high values going from purple to 
red. Linestyles are as in Figure [3] 



cases, with a slight reduction in the average from 1.08733 for 
the isothermal a priori to 1.05709 for the adiabatic a priori, 
so the choice of a priori does not affect the quality of fit 
to the spectrum. In regions of low information content, the 
retrieved temperature profile will be affected by the a priori 
and therefore it is important to perform retrievals from a 
range of different a priori starting points to ensure that this 
behaviour is fully understood. 

Similar tests for 100 spectra have been performed in 
transmission geometry; these tests are complicated by the 
fact that there is insufficient sensitivity to temperature in 
transmission spectra to disentangle the effects from those 
due to gaseous absorption, so in order to perform a retrieval 
one must assume a temperature structure. We test the ef- 
fect of incorrect assumptions about the T-P profile on the 
retrievability of the gas VMRs. First, we assume the same 
constant temperature profile as the a priori for the eclipse 
case. It is reasonable to assume a constant temperature that 
is within the range of temperatures in the 'real' profile, as it 
is possible to compute an equilibrium temperature (Equa- 
tion [3]) for an exoplanet based on measurement of its orbital 
parameters: 



(1-^) 



4a2 



(3) 



where -k indicates stellar radius R and temperature T, 
A is the Bond albedo of the planet and a is the orbital semi- 
major axis. This is often a poor approximation for giant 
planet temperatures in our own solar system, as they also 
have large internal heat sources, but is the best first-order 
approximation possible for extrasolar planets. 

Secondly, we perform the retrieval with the tempera- 
ture profile from the forward model that we used to create 
the noisy synthetic. Since we could use a temperature con- 
straint from eclipse measurements as an input for a trans- 
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Figure 6. Comparison of temperature retrievals from three dif- 
ferent a priori starting points (top panel). The mean and stan- 
dard deviation over 100 retrievals are shown for each case; the 
joined crosses indicate mean values and the shaded envelope the 
standard deviation. The a priori temperatures are shown on the 
bottom. The input temperature profile is the overplotted black 
solid line in both panels. 

mission retrieval, which with perfect heat redistribution in 
the atmosphere would be the same as the temperature at the 
terminator, it is conceivable we could approach this level of 
knowledge about the terminator temperature; the level of 
prior knowledge about temperature structure for the major- 
ity of cases is likely to fall somewhere between these extrema. 

The results for H2O, CO2 and CH4 in these two cases 
are shown in Figures [71 and [H] results for CO and NH3 are 
as poor as the results obtained from eclipse spectra. The 
results presented for the no-inversion case have assumed a 
noise level consistent with combining measurements from 
30 transits to increase the signal to noise, whilst noise for a 
single transit has been added to the inversion case. Adding 
together 30 transits improves the retrieval accuracy for cases 
where the VMRs are low, provided the assumed temperature 
profile is close to reality, but no improvement is observed for 
the constant temperature profile. However, the difference in 
/t^ between the correct and incorrect temperature profile 
assumptions is much greater when the measurement error is 
smaller, so improving signal to noise increases the likelihood 
of being able to use the to distinguish between different 
temperature scenarios. This can be seen in Figure O the fit 
is poor in both cases (30 observations and 1 observation) 
where the temperature is assumed to be isothermal, but the 
difference in goodness of fit is much easier to discern in the 
left hand panels, where the signal to noise is higher. 

As can be seen from these plots, an assumption of con- 
stant temperature results in a very poor retrieval of H2O 
VMR and significantly reduces the retrieval accuracy of CO2 
and CH4. However, the average /'^ for the constant tem- 
perature case is 30, whereas it is close to 1 when the cor- 
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Figure 7. Retrievals of H2O, CO2 and CH4 from transmission spectra for a hot Jupiter with no temperature inversion. In the upper 
panels, a constant temperature of 1500 K is assumed; in the lower panels, the input temperature profile is assumed, colours are as 
Figure [3l 
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Figure 8. As Figure [71 but for a hot Jupiter with a temperature inversion. 



rect temperature profile is used. Although the temperature 
structure cannot be directly retrieved from a transmission 
spectrum, it may therefore still be possible to place some 
constraint on it by searching for a minimum within a 
range of different assumed temperature profiles. This could 
be done by repeating the retrieval over a range of differ- 



ent temperature assumptions and choosing the solution with 
minimum x^. 

Even if we have no temperature information from a sec- 
ondary eclipse measurement, we can still do better than the 
isothermal assumption by calculating the expected equilib- 
rium temperature Teq of the planet and modelling the at- 
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Figure 9. Examples of fits to noisy synthetic hot Jupiter transmission spectra. The 1-cr error bar is shown in dark gray shading, and the 
2-cr error in fight gray. The joined crosses in the top panel show the noisy synthetic, and the solid black line shows the fitted spectrum. 
The left-hand panels are the hot Jupiter spectra without a temperature inversion (30 transits) and the right show the Jupiter with a 
temperature inversion (single observation). The red shaded areas indicate the regions where the spectral fit is poor when an isothermal 
temperature approximation is used. 



mosphere as a single slab. This slab of atmosphere radiates 
equally out to space and down into the lower atmosphere at 
the stratospheric temperature Tstrat- If we equate the total 
radiation absorbed and emitted by the slab, we can relate 
these two temperatures using Equation [4] 



(4) 



where e is the emissivity of the slab, a is the Stefan- 
Boltzmann constant and the factor of 2 arises from the fact 
that the slab receives radiation from space only but emits 
radiation equally in both directions. Rearranged, this gives 



the relation that Tec 



We can then calculate the temperature at the bottom of 
the troposphere/surface of the planet by calculating the dry 
adiabatic lapse rate of the atmosphere, given by Equation [S) 
g is the acceleration due to gravity and Cp the specific heat 
capacity. This is to some extent arbitrary as we do not know 
a priori the pressure level at which this occurs, but for a hot 
gaseous planet we assume that this occurs at approximately 
the 1 bar level; the retrieval is insensitive to temperature 
variations below this level (see Figure [J). The top of the 
tropopause we assume to occur at around 0.1 bar, based on 
observations of solar system giant planets. We can assume 
that pressure and altitude are related by the hydrostatic 



equation, [G] 



r = 



dT _ 



9_ 



(5) 



(6) 



H is the atmospheric scale height, which is given by 
H — kT/mg where k is the Boltzmann constant, T is tem- 
perature and m is the molecular mass of the atmosphere. 
If we know the bulk atmospheric composition and assume 
for the purpose of calculating the scale height that T is 
the stratospheric temperature, we can combine equations [5] 
and [6] to calculate the temperature at the base of the tropo- 
sphere. 



^ trop 



Tstr: 



mg \p2 



(7) 



Using this equation, we can derive a simple temperature 
profile based on the planet's equilibrium temperature. 

We compare the results from all three temperature pro- 
file parameterisations for a hot Jupiter without a tempera- 
ture inversion in Figure [10] Even though the shape of the 
original temperature profile is better reproduced by the sin- 
gle slab atmosphere model, it does not significantly improve 
the retrieval result over the isothermal assumption case. The 
H2O retrieval is improved when the estimated temperature 
structure is more realistic, but the CO2 and CH4 retrievals 
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A comparison between transmission retrievals for the hot Jupiter case where different temperature profiles (shown in the 
assumed. The joined crosses are the mean values and shaded areas represent the standard deviation over 100 retrievals. 



are not as good as those obtained with the isothermal as- 
sumption. The /n is also higher for the single slab model 
than it is when the correct temperature profile is used. 

The interpretation of transmission measurements with- 
out a good temperature constraint from eclipse is therefore 
something to be undertaken with caution. Even if eclipse 
results are available and the dayside temperature is con- 
strained, it is important to bear in mind that hot, short- 
period planets are often tidally locked, so we might ex- 
pect some variation in temperature between the limb as 
seen in transmission and the star-facing disc as seen in 
eclipse; any assumptions made about heat redistribution will 
also affect the retrieval of gaseous abundances in transmis- 
sion geometry. Some information on this score can be ob- 
tained using phase curves/transit mapping, as in the anal- 
yses of HP 189733b performed bv lMaieau et al and 
Ide Wit et all (l2012l ). which provides information about the 
brightness temperature as a function of local time on the 
planet. The limitations of optimal estimation mean that it 
is not possible to use something like NEMESIS to indepen- 
dently constrain temperature and gaseous abundances from 
a primary transit spectrum, but the difference in /n be- 
tween the different scenarios indicates that there is some 
information available in the spectrum about the tempera- 
ture. This can be exploited by performing multiple retrievals 
with different temperature profiles and comparing the i as 
suggested above. 

We also find that there is degeneracy between the as- 
sumed planetary radius and the temperature in retrievals 
of transmission spectra. In Figure [TTl we show /n for a 
range of model fits where the planetary radius and temper- 
ature are perturbed from the known input value. It is clear 
that for a range of radii and temperatures, ±150 km and 
± 60 K respectively, a reasonable fit to the noisy synthetic 
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Figure 11. A map of x^/^ when the radius and temperature of 
the best-fit model to a noisy synthetic are perturbed from the 
correct value. There are clearly several possible solutions with 
X^/n of less than 2, indicating that there is degeneracy between 
planetary radius and temperature. 

is obtained. This degeneracy arises from the fact that in- 
creasing the planet radius and increasing the temperature 
have the same effect on a transmission spectrum, as both 
make the planet appear to be larger with an increased scale 
height. Therefore, a simpler way to account for a lack of in- 
formation about temperature structure might be to adjust 
the modelled radius until /n is minimized; this would al- 
low us to obtain accurate information about the atmospheric 
composition even if we do not know the temperature. 

4.2 Warm Jupiter 

We now consider the case of a somewhat cooler Jupiter- 
size planet. The temperature structure is similar to that of 
the hot Jupiter without a temperature inversion, but the 
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Figure 12. Retrieval results for a warm Jupiter with a single eclipse; colours/linestyles as Figure [3] 
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Figure 13. Comparison of retrieval results for a warm Jupiter assuming errors for a single eclipse (black/grey) and 30 observations 
combined (turquoise/blue). Joined crosses are the mean over 100 retrievals, shaded areas are standard deviations. The solid black lines 
represent the input temperature and the 1:1 correspondence line in the temperature and gas plots respectively. The increase in SNR for 
30 transits much improves the quality of the retrieval. 



temperature at 1 bar is 1000 K instead of 2500 K. This 
reduces the planet/star contrast and the signal to noise in 
eclipse, and in Figure [12] we see that the retrieval for a single 
eclipse is poor. 

The signal to noise could be improved in this case if sev- 
eral observations were added together. We repeat the test 
with the added noise reduced by a factor of a/30, which is 



the equivalent of observing 30 transits and combining the re- 
sults assuming that all error is random rather than system- 
atic; as expected, the quality of the retrieval is significantly 
improved (Figure [13]). The spectra for both cases are shown 
in Figure [14] Based on the input temperature profile this 
planet would have an equilibrium temperature of around 
600 K, which corresponds to an orbital radius of 0.3 AU for 
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Figure 14. Fits to noisy synthetic warm Jupiter eclipse spectra. The residuals are shown in the lower panels. The 1-cr error bar is shown 
in dark gray shading, and the 2-cr error in light gray. The joined crosses in the top panels show the noisy synthetic, and the solid black 
lines show the fitted spectra. 




Figure 15. Retrievals of H2O, CO2 and CH4 from single transmission spectra for a warm Jupiter with no temperature inversion. In the 
upper panels, a constant temperature of 600 K is assumed; in the lower panels, the input temperature profile is assumed. Colours are as 
Figure [3] 



a sun-like star; it would therefore take around 60 days to 
complete one orbit. If at least 30 transits are required be- 
fore sufficient signal to noise is achieved, then such a planet 
would require observation over EChO's proposed five-year 
lifetime. This means that retrievals of similar planetary at- 
mospheres are a possibility, but the periods of planets like 
this would need to be well characterised prior to the start 



of the mission, perhaps using ground-based telescopes, to 
ensure that the efficiency of observation is maximised. It is 
wo rth noting that vary ing stellar activity may, as mentioned 
bv lTinetti et al.1 (I2OIOI ). complicate the process of averaging 
over several observations. 

The quality of transmission observations is also reduced 
for a cooler planet (Figures [15] and [16]); the magnitude of 
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the transmission signal depends on the scale height of the 
atmosphere, which is proportional to temperature, so the 
warm Jupiter has a compressed atmosphere relative to the 
hot Jupiter. Example spectra for two different temperature 
assumptions are shown in Figure [TTl along with those for 
the hot Neptune case. As with the hot Jupiter, a good prior 
constraint on the temperature is necessary for a reasonable 
retrieval of gaseous abundances, so it will probably be nec- 
essary to perform multiple retrievals with different tempera- 
ture profiles and search from the minimum x^, as suggested 
for the hot Jupiter. A good temperature constraint is harder 
to obtain from eclipse measurements for a cooler planet, and 
the difference in /n between good and poor solutions is 
smaller for spectra with larger error bars, so in both transit 
geometries cooler planets present an observational challenge. 

4.3 Hot Neptune 

Finally, we consider a planet approximately the size of Nep- 
tune orbiting an M dwarf. Even though this planet is smaller 
and cooler (1500 K in the deep atmosphere) than the hot 
Jupiter we consider, because it is orbiting a smaller, cooler 
star the contrast ratios are better (Figure [TH]). We can suc- 
cessfully retrieve the temperature profile and VMRs of H2O, 
CO2 and CH4 assuming noise levels for a single transit (Fig- 
ure [20}. We can also retrieve VMRs of CO greater than 10 
ppmv. 

The hot Neptune has an extended H2-He atmosphere 
and is orbiting a relatively small star, so the SNR in trans- 
mission is very high (Figure [19]). As for the hot Jupiter 
cases, when the assumed temperature profile is correct the 
retrievals of H2O, CO2 and CH4 VMRs are extremely accu- 
rate, but when an isothermal temperature profile is used the 
retrieved value is incorrect and the fit is poor (Figures [2T1 
and 1191 ). As with the hot Jupiter, the small errors on the 
spectrum mean that poor assumptions about the tempera- 
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Figure 17. Examples of fits to a noisy synthetic warm Jupiter 
transmission spectrum, (as Figure [Oj. 
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Figure 18. Fits to a noisy synthetic hot Neptune ecHpse spec- 
trum, as Figure [M] 

ture profile are associated with high x^, so it should be pos- 
sible to achieve a good retrieval of gas properties by varying 
the assumed temperature and repeating the retrieval. 



5 DISCUSSION 

For the simple model planets described above, our results in- 
dicate that it will be possible for the EChO space telescope, 
and any other telescopes with a similar spectral range and 
an SNR/spectral resolution as good as EChO's, to reliably 
retrieve temperature as a function of pressure and VMRs for 
H2O, CO2 and CH4, with upper limits for CO and NH3. We 
now consider the validity of some of the simple assumptions 
made here, and discuss how changes in these might affect 
our conclusions. 

5.1 Planetary mass and radius 

We have assumed that we have perfect knowledge of the 
planetary mass and radius. The methods used to calcu- 
late these parameters for transiting exoplanets mean that 
their precision depends on the precision to which we know 
the mass and radius of the parent star. Errors quoted in 
iDemorv et al.l (|20Q9l ) for Proxima Centauri (GJ 551) are ap- 
proximately 5 %, therefore we may expect at least a 5 % er- 
ror in measurements of planetary mass and radius. We also 
assume the measured transit radius is equivalent to the ra- 
dius at the bottom of our atmosphere, i.e. at 10 bar pressure, 
which may not be the case in reality as the atmosphere is 
likely to be opaque to higher altitudes in some wavelengths. 
The absolute radius derived from the transit depth averaged 
over all visible wavelengths of our model hot Jupiter is 2% 
greater than the specified radius at 10 bar. 

The worst case scenario is that the derived mass is too 
large/small when the derived radius is too small/large, as 
this produces the largest error in the calculated gravitational 
acceleration. The scale height and dry adiabatic lapse rate, 
and therefore the temperature profile, depend on the grav- 
itational acceleration, so large variations in g can make a 
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Figure 19. Examples of fits to a noisy synthetic hot Neptune 
transmission spectrum, (as Figure [9]). 

significant difference to spectra. We test the effect of this 
by comparing retrievals assuming 1) the actual input radius 
and mass, 2) 95% of the mass and 105% of the radius and 
3) 105% of the mass and 95% of the radius. The retrieved 
temperature profiles in eclipse are compared below, along 
with the /n values. 

These extreme cases of inaccuracy in the planetary mass 
and radius make a large difference to the retrieved temper- 
ature, CO2 and CH4 values. A too-low gravity has a more 
significant effect that a too-high gravity. Using the wrong 
values for mass and radius increases the parameter to 

a value above 1, which indicates that the fit to the noisy 
synthetic spectra becomes poor; the retrieved temperature 
profiles are also less smooth, particularly for the case where 
the gravity is low. This means that it should be possible 
to reject a result where the wrong mass and radius are as- 
sumed for a real planet, and indeed it may be possible to 
constrain the correct range of g by minimising /n over 
different values of mass and radius. 

We also test the impact of errors in the mass and radius 
on results in transmission geometry. The same deviations 
as with the eclipse case mean that NEMESIS is unable to 
produce a fit to the spectrum at all, resulting in x^ /n values 
of several thousands. As with previous results of a similar 
nature, the extremely high values of x? /n indicate that such 
a retrieval would not be accepted, and so even if errors on 
the mass and radius of the planet are present it should still 
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Figure 20. Eclipse retrieval results for a hot Neptune orbiting an M star; colours/linestyles as Figure [S] 
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Figure 21. Transmission retrieval results for a hot Neptune orbiting an M dwarf (as Figure [7}. 



be possible to obtain sensible constraints through trial and 
error of different planetary properties. However, as shown 
in Figure [TT] there is significant degeneracy between radius 
and temperature in primary transit, and we would expect 
degeneracy between g and temperature too as both affect 
the atmospheric scale height {H = kT/mg), so constraining 
g in this way without any constraint on temperature would 
be difficult for primary transits. 

Even if errors in the planetary mass and radius are as 



great as 5% when these are derived from stellar mass and 
radius, it is possible to place more precise constraints on 
the planetary surface gravity by calculati ng; it directly from 
trans it and radial velocity observables (jSouthworth et all 
I2OO7I ) ■ A large source of error in planetary mass and radius 
is a lack of constraint on stellar mass and radius, whereas 
this alternative method does not require knowledge of these 
parameters. The large effect of uncertainty in g on the re- 
trievals discussed here means that the method of Southworth 
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Figure 22. A comparison between retrievals with the correct planet mass and radius (black/grey) and retrievals with a low mass/high 
radius (blue/turquoise) and a high mass/low radius (red/fuschia). The joined crosses are the mean values and shaded areas represent 
the standard deviation over 100 retrievals. The yi^ jn is also shown. 



et al. is likely to be very important for calculating the grav- 
itational acceleration of new planet candidates. 



5.2 Systematic noise 

As mentioned in Section [331 the analysis so far has assumed 
photon-limited noise, which is a perfect case; in reality, we 
expect some systematic contribution to the noise from the 
instrument and from sources such as stellar variability. It is 
a useful exercise to investigate the effect of multiplying the 
photon noise by a range of factors, and determine at what 
point the signal is lost for each of the cases above. The re- 
sults of these tests are shown in Figures [23] — [25] for the 
hot Jupiter and hot Neptune cases, where the colours cor- 
respond to the increase in the noise level. Noise levels up 
to 40 X the photon noise were tested in order to investigate 
the point at which the retrieval completely fails, but more 
likely values are in the 2 — 4x range. It can be seen that in 
eclipse, a factor 4 increase (very dark purple) from the pho- 
ton noise limit does not substantially affect the quality of the 
temperature retrieval, and a factor 2 increase (black) does 
not affect the gas volume mixing ratio retrieval for higher 
abundances. In transmission, the effect for a factor 2 in- 
crease is not significant but a factor 4 increase does start to 
affect the retrieval accuracy. If however a factor 3 increase 
represents a worst-case scenario (see Section 13. 3|) , this indi- 
cates that the inclusion of a reasonable level of systematic 
noise does not substantially change the conclusions we ar- 
rive at. For a factor 2 increase in the noise level for the hot 
Jupiter, we can only detect H2O down to an abundance of 1 
ppmv, CO2 down to 5 ppmv and CH4 down to 0.5/5 ppmv 
(eclipse/transmission); for the hot Neptune, the detection 
limit for a factor 2 increase is unchanged, but for a factor 4 
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Figure 23. Effect of increasing the eclipse spectrum noise on 
the temperature retrieval for a hot Jupiter with a temperature 
inversion, orbiting a G2 star (upper panels) and a hot Neptune 
orbiting an M5 star (lower panels). Colours correspond to the 
noise level relative to the photon-limited case. 



it increases to 5 ppmv in transmission and eclipse for H2O 
and CO2, and 5 ppmv for CH4 in eclipse. At around 30 x 
the photon noise there is not enough signal remaining for 
the retrieval to move away from the a priori. 

We have not shown the effect of increasing the noise on 
the warm Jupiter spectra because we are already consider- 
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Figure 24. Effect of increasing the eclipse spectrum noise on the gas abundance retrieval for a hot Jupiter with a temperature inversion, 
orbiting a G2 star (upper panels) and a hot Neptune orbiting an M5 star (lower panels). 




ing adding together multiple transits to gain sufficient signal 
to noise; the effect of systematic noise would simply be to 
increase the number of transits required to achieve a given 
retrieval accuracy, although in this case a small noise in- 
crease may result in the mission lifetime being insufficient for 
the number of observations required. However, the averaging 
process introduces another source of systematic noise, that 
of changing stellar activity. The majority of stars have spots 
on their surfaces, which tend to be cooler than the rest of 



the disc and also have different spectral signatures. The frac- 
tional spot coverage is unlikely to change on the timescale of 
a single transit, but between visits of around 60 days apart, 
as for the warm Jupiter, we may expect it to vary. Star spots 
can introduce error into planetary transmission spectra be- 
cause the part of the stellar disc o cculted by the pla net may 
not be representative of the whole. IPont et al.1 (|2Q12l) discuss 
the effects of unocculted star spots on transmission spectra. 
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Figure 26. Warm Jupiter transmission results over 30 transits 
where the star is unspotted (black points) and where spot-like 
effects have been added to the synthetic spectrum (red points). 
Both retrievals have assumed that the star is unspotted. It can 
be seen that star spots are an extra source of error, but retrievals 
are still possible. 



and show that the magnitude of the effect is most significant 
at optical wavelengths. 

Assuming that the ratio of the spot spectrum to the 
unspotted solar spectrum can be approximated as a ratio of 
Planck functions at 4800 and 5800 K respectively, we use the 
formulae given in IPont et al.1 (l2012l ) to calculate the adjusted 
planet/star radius ratio for a range of spot flux dimmings 
at 0.6 i^m. We a ssume, based on the magnitudes given in 
IPont et al.1 (l2012l ). that the la flux dimming is 0.25 % at 
0.6 i^m, and we generate 30 spectra with different amounts 
of star spot dimming, plus photon noise, for each test atmo- 
spheric composition. We then average over these 30 spec- 
tra and perform a retrieval on ten examples, enabling us to 
quantify the increase in error due to stellar variability. The 
results of this test are shown in Figure [26l It can be seen that 
the error on the retrieved H2O and CH4 VMRs increases by 
around a factor 3, whereas the error on the CO2 VMR could 
be as much as an order of magnitude. However, no attempt 
has been made to account for the spot effect in the retrieval, 
so this could conceivably be improved upon if the stellar flux 
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Figure 27. Variation in retrieved temperature if H2O 
(black/grey), CO2 (blue/turquoise) and CH4 (red/fuschia) are 
not included in the retrieval model. The joined crosses are the 
mean values and shaded areas represent the standard deviation 
over 100 retrievals. The x^/n is also shown for each case. 

is monitored and the inf e rred s pot coverage corrected for, as 
described in lPont et all (l2012f ). 

If only one or two hot Neptunes around M dwarfs are 
observed by EChO, it is likely that these would also be ob- 
served several times. We perform a similar test for the sen- 
sitivity to star spots in this case, with a spot temperature 
of 2000 K compared with a stellar temperature of 3000 K, 
and find that the effect is negligible due to the much larger 
signal in transmission than that of the warm Jupiter. 

We do not consider occulted spots in this analysis be- 
cause these should have observable signa tures in the visible 
wavelength lightcurves (|Pont et al.ll20lj V so are somewhat 
easier to correct for. We also do not consider the effect of 
star spots on eclipse measurements because we do not ex- 
pect significant spot evolution during an eclipse, and as the 
full stellar disc is visible throughout the spot distribution 
does not affect the measurement. 



5.3 Missing absorbers 

We test the effect of removing gaseous absorbers from the 
retrieval model that were included in the synthetic spectrum 
calculation. We perform this test for H2O, CO2 and CH4 as 
these are the gases that produce the largest effect on the 
spectra, as evidenced by the reliability of the retrievals for 
these gases. If any one of these gases is ignored, the tem- 
perature retrieval becomes less accurate, with H2O having 
the most significant effect. The /n value is increased from 
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Figure 28. Top: 100 CH4 retrievals assuming constant mixing ra- 
tio. Colours correspond to jn as previously. The input values 
are the VMRs at 0.1 bar. Bottom: 6 examples of a continuous CH4 
retrieval. The enclosed shaded areas represent the error envelope 
of each retrieval. The dotted profiles of the same colours are the 
corresponding model inputs. The black dashed line is the a •pri- 
ori profile. The hatched areas indicate the atmospheric pressure 
regions for which the CH4 VMR can be retrieved in the majority 
of cases. 



the values when all the gases are included, but the correct 
shape of the temperature profile is still retrieved. As gas ab- 
sorption line data at high temperatures is often incomplete, 
it is necessary to be aware that missing lines can impact the 
reliability of temperature retrieval. 



5.4 Altitude variation of absorbers 

We also test the effect of including m ore realistic p rofiles 
of gas VMR as a function of altitude. iMoses et al.l (|201lh 
use a ID photo- and thermo-chemical model to calculate the 
expected abundances of trace absorbers on HD 189733b and 
HD 209458bb. They find that CH4 and NH3 are expected 
to be photochemically removed at high altitudes. We test 
the effect on our results of including a more realistic vertical 
CH4 profile, based on the prediction of Moses et al. (201 J). 

We first test the effect of including a more complex CH4 
profile in the model but performing the retrieval in the same 
way as before, i.e. only retrieving a single multiplying factor 
on a constant VMR. Despite this crude assumption, we can 
accurately retrieve the abundance of CH4 at 0.1 bar (Fig- 
ure [28l upper panel), which is within the altitude range of 
maximum sensitivity. Ignoring the detail of the CH4 vertical 
profile has very little effect on the retrieval of other quanti- 
ties, except for increased error on the temperature close to 
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Figure 29. 100 temperature retrievals when a more realistic CH4 
vertical profile is included in the model atmosphere. In the upper 
panel is the temperature retrieved when a constant CH4 profile 
is still assumed in the retrieval; in the lower is the temperature 
retrieved when CH4 VMR is retrieved as a function of pressure. 
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Figure 30. As Figure [28] but for transmission. 
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Quantity 


Hot Jupiter 


Warm Jupiter 


Hot Neptune 


Temperature 


200 K (1 bar— 1 mbar) 


250 K (0.5 bar— 5 mbar) 


100 K (1 bar— 1 mbar) 


H2O 


2x (>0.5 ppmv) 


4x (>50 ppmv) 


2x (>0.5 ppmv) 


CO2 


3x (>10 ppmv), 6x (>0.1 ppmv) 


5x (>50 ppmv) 


2x (>0.1 ppmv) 


CO 


5x (>100 ppmv) 


Not detectable 


3x (>10 ppmv) 


CH4 


2x (>1 ppmv), 3x (>0.1 ppmv) 


4x (>50 ppmv) 


1.5x (>0.1 ppmv) 


NH3 


5x (>100 ppmv) 


Not detectable 


3x (>100 ppmv) 



Table 3. Information available from a single observation for each of the test planets in eclipse. 
These precisions are only valid for certain pressures (for temperature) and concentrations (for 
gases), which ranges are given in parentheses for each case. Precisions are quoted as zbK for 
temperature. As the gas abundances vary over several orders of magnitude, the precisions are 
given as multiples of the abundance; for example, the H2O VMR can be retrieved to within a 
factor 2 of the true value for a hot Jupiter, provided there is more than 0.5 ppmv present. The 
values quoted take into account the errors due to retrieval degeneracy, which is the dominant 
source of error, by calculating the precision based on the 2cr deviation from the true value over 
all retrievals. 



Quantity 


Hot Jupiter 


Warm Jupiter 


Hot Neptune 


H2O 


1.5x (>0.5 ppmv) 


1.5x (>10 ppmv) 


1.5x (>0.1 ppmv) 


CO2 


1.5x (>0.5 ppmv) 


2x (>50 ppmv) 


1.5x (>0.1 ppmv) 


CO 


Not detectable 


Not detectable 


1.5x (>10 ppmv) 


CH4 


1.5x (>0.5 ppmv) 


2x (>10 ppmv) 


1.5x (>0.1 ppmv) 


NH3 


3x (>100 ppmv) 


Not detectable 


1.5x (>5 ppmv) 



Table 4. Information available from a single observation for each of the test planets in transmis- 
sion. Perfect knowledge of the temperature profile has been assumed. Otherwise everything is as 
Table [31 



the 10 bar level, which is reassuring as it implies little degen- 
eracy between the sort of altitude variation in gas abundance 
that we would expect and changes in temperature structure 
(Figure [291 

upper panel). 

We also investigate the possibility of obtaining informa- 
tion about the vertical profile of CH4 from eclipse measure- 
ments by performing a continuous retrieval of CH4 VMR 
as a function of pressure. Whilst this removes the error in 
retrieved temperature in the deep atmosphere (Figure [29l 
lower panel) , it is apparent that not much information about 
the CH4 profile is available from the spectrum (Figure I28[ 
lower panel). The continuous retrieval only matches the in- 
put at two pressure levels, and everywhere else relaxes back 
to the a priori solution, so we do not obtain much more use- 
ful information by increasing the complexity of the retrieval 
even if absorber abundances are variable with altitude, /n 
is also slightly higher on average (1.1 rather than 1.0) for the 
continuous retrieval versus the constant VMR case because 
retrieving a continuous profile of CH4 reduces the number 
of degrees of freedom; a continuous retrieval has more in- 
dependently varying parameters than a multiplying factor. 
This indicates that the improvement in the fit to the spec- 
trum, if any, is not sufficient to justify the increased number 
of retrieved parameters. Therefore, the simpler approach of 
assuming a vertically constant VMR for all gases is preferred 
for eclipse retrievals. 

We perform the same test for the transmission case 
(Figure [30]) . The results are similar to those for eclipse, with 
a limited amount of vertical information available, but the 
altitudes of sensitivity are slightly higher in transmission 
compared with eclipse; a combination of transmission and 
eclipse measurements could provide enough information to 



constrain variation with altitude, provided the temperature 
variation between the dayside and terminator regions is not 
significant. In both cases, the retrieval ceases to be accurate 
for tropospheric VMRs lower than 10~^, so for values below 
this an upper limit only would be achievable. 



6 CONCLUSIONS 

We find that for a hot Jupiter orbiting a sun-like star and 
a hot Neptune orbiting a small M-dwarf a single eclipse ob- 
servation by an EChO-like telescope is sufficient to provide 
good constraints on temperature structure and volume mix- 
ing ratios of H2O, CO2 and CH4. Some constraint is also 
possible for CO and NH3 VMRs, although when these val- 
ues are low an upper limit only is achievable. The lower 
sensitivity to temperature structure in transmission obser- 
vations means that independently constraining the temper- 
ature and gas VMRs from transmission spectra is not pos- 
sible, but nonetheless errors in the assumed temperature 
structure do have a significant effect on the retrieved gas 
VMRs. The constraints on the gas VMRs are good provided 
the temperature structure has already been accurately re- 
trieved from eclipse observations, so transmission and eclipse 
spectra combined could provide information about poten- 
tial changes in atmospheric composition between the day- 
side and terminator. This approach of course relies on the 
assumption that the temperature profile at the terminator 
is the same as, or at least similar to, the dayside averaged 
temperature profile that is retrieved from the eclipse spec- 
tra; whether or not this is the case depends heavily on the 
efficiency of heat transport from the day to the night side of 
the planet, which is not known a priori. A combination of 
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transit mapping (|de Wit et al.ll2012h . secondary eclipse and 
transmission measurements from a single orbit would there- 
fore provide a wealth of information about hot Jupiters and 
Neptunes. These findings are summarised in Tables [3] and [H 
We find that there is little dependence on the cho- 
sen a priori temperature profile for eclipse retrievals, pro- 
vided there is information contained in the spectrum; this 
demonstrates the reliability of the retrieval method used. 
We also find that retrievals are not very sensitive to verti- 
cal variation in the abundance of CH4, a gas expected to 
vary as a function of pressure on planets like HD 189733b 
(|Moses et al.ll201ll ). justifying the common assumptions of 
constant VMRs for trace absorbers. The omission from the 
retrieval model of absorbers that are present in the origi- 
nal model atmospheres does however reduce the accuracy 
of the retrieved temperature profile, suggesting the need for 
caution in interpreting results where absorption line data 
is known to be incomp lete. Theoretical programmes such 



as the Exomol project (iTennvson Yurchenko 



new experimental results (e.g. iHargreaves et al 



20111 



2012) and 



increase the robustness and reliability of retrieval results by 
improving the available line data. 

Errors in the measured radius and mass of exoplanets 
could have an adverse effect on retrieval quality. A 5% error 
in the mass and a 5% error in the opposite sense in the 
radius result in the retrieved temperature profile in eclipse 
being wrong by as much as 500 K, which in turn impacts 
the retrieval of trace gas abundances. Similar errors in mass 
and radius make it impossible to produce a reasonable fit to 
transmission spectra. The high values for these cases 

indicate that the retrievals are not reliable, so it may still be 
possible to fit the spectra well by re-determining the mass 
and/or radius through minimising /n by trial and error. 

Whilst the level of systematic noise expected for EChO 
is not yet known, we investigate the effect of increasing the 
noise from the photon-limited case and find that good tem- 
perature retrievals are still possible for the hot Jupiter and 
Neptune cases even with a factor of 4 increase, and good gas 
VMR retrievals can be obtai ned for a fact o r of 2 increase. 
Based on the conclusions of iGibson et al.1 (l201lh . we con- 
sider a factor of 3 increase over the photon-noise limit to be 
a worst case scenario, so do not expect systematic noise to 
greatly alter our conclusions. 
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